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ABSTRACT 



Context. We report the results of our XMM-Newton monitoring of SN1987A. The ongoing propagation of the supernova blast 
wave through the inner circumstellar ring caused a drastic increase in X-ray luminosity during the last years, enabling detailed high 
resolution X-ray spectroscopy with the Reflection Grating Spectrometer. 

Aims. The observations can be used to follow the detailed evolution of the arising supernova remnant. 

Methods. The fluxes and broadening of the numerous emission lines seen in the dispersed spectra provide information on the evolution 
of the X-ray emitting plasma and its dynamics. These were analyzed in combination with the EPIC-pn spectra, which allow a precise 
determination of the higher temperature plasma. We modeled individual emission lines and fitted plasma emission models. 
Results. Especially from the observations between 2003 and 2007 we can see a significant evolution of the plasma parameters and a 
deceleration of the radial velocity of the lower temperature plasma regions. We found an indication (3cr-level) of an iron K feature in 
the co-added EPIC-pn spectra. 

Conclusions. The comparison with Chandra grating observations in 2004 yields a clear temporal coherence of the spectral evolution 
and the sudden deceleration of the expansion velocity seen in X-ray images ~6100 days after the explosion. 

Key words. ISM: supernova remnants 
- Shock waves 
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1. Introduction 

The circumstellar ring system around SN 1987 A was ejected by 
the progenitor star approximately 20000 years before the su- 
pernova explosion. About 10 years after the explosion the blast 
wave started to propagate through the inner ring, causing com- 
pression, heating and ionization of its matter At that time bright 
unresolved regions, so called hot spots, ap peared all around the 
ring in HST images (iLawrence et al ] l2000l) and the rather linear 
increa se in the soft X-ray band, as observed with RO SAT since 
1992 (iBeuermann et al.l |l994- Hasi nger et all 1 19961) . was fol- 
lowed by an exponential brightening as monitored with XMM- 
Newton, Chandra, Suzaku an d Swift. Recently, some flattening 
of the flux increase is seen (ISturm et al.ll2009l and references 
therein for the individual fluxes). 

The X-ray spectra are interpreted as thermal emission com- 
posed of a lower temperature component (~0.5 keV) and a 
higher temperature component (~2.5 keV). Deep Chandra grat- 
ing observations allowed the measuring of the bulk gas veloc- 
ity due to spectral line deformation. Surprisingly these values 
were lower than the velocities expected from the plasma temper- 
atures. This suggests the contribution of reflected shocks, addi- 
tional to the 'normal' for ward shock (Zhekov et al. 2005, 2006, 
12009: Dewev et al.ll2008h . The expansion velocity derived from 
Chandra X-ray images is even higher. About 6100 days after the 
explosi on a sharp decelera tion to 1600-2000 km s ' was ob- 
served ("Racusin et al."2009^. 

Shocks transmitted into denser regions of the ring have a 
slower shock wave velocity and therefore can be responsible for 
the low temperature component. This interpretation is supported 
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by th e morphology seen in optical and X-ray images dMcCravl 
l2007 l). as well as the similar evolution of the soft X-ray flux (0.5- 
2.0 keV) and the evoluti on of highly ionized opti cal emission 
lines from the hot spots (iGroningsson et al1l2006h . This emis- 
sion might be caused by even slower radiative shocks. 

R enewed radio emis sion of SN 1987 A was detected in July 
1990 (iTurtle et all 19901) . A continuously rising flux and increas- 
ing source radius has been obs erved by the ATCA since then 
(iGaensler etalJ 12001 iNg etalj r2008). The increase of the ra- 
dio light curve match es the evolu tion of the hard X-ray flux (2- 
10 keV) quite closelv (iPark et al.l 2005: Aschenbach 2007), thus 
the synchrotron radio emission may originate in the hot thermal 
plasma between the forward and reverse shock and a fraction of 
the hard X-ray flux may also have a non-thermal origin. Broad 
Ha and Lya lines suggest the presence of the reverse shocks 
(Michael et al. 2003; Heng et al. 2006). 

This study reports on the yearly XMM-Newton monitoring 
observations of SN 1987 A between January 2007 and January 
2009 together with one prior observation from May 2003, mak- 
ing it possible to reveal the detailed plasma evolution. We mea- 
sured light curves and widths of individual emission lines and 
analyzed the spectral plasma evolution by fitting thermal plasma 
emission models. 



2. Observations and Data Reduction 

Since January 2007 we have per formed a yearly mo nitoring of 
SN 1987 A using XMM-Newton ('Jan sen et al1l200ll) . Including 
an earlier observation from May 2003, four observations with an 
exposure each exceeding 100 ks were accomplished. A summary 
is listed in Table [1] Earlier observations and the May 2003 data 
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Table 1. XMM-Newton observations of SN 1987 A and data selection. 



Date 


Time 
(UT) 


Satellite 
Revolution 


Instrument 


Read-out 
Mode 


Filter 


Net Exp 
[s] 


Counts*"' 


2003 May 10-11 


16:41-23:37 
11:42-23:38 
10:47-23:38 


0626 


EPIC-pn 
RGSl 
RGS2 


SW, 6ms 
Spectro 
Spectro 


medium 


34983 
112164 
111671 


13712 
3061 
3635 


2007 Jan 17-19 


19:28-01:18 
18:23-01:19 
18:23-01:18 


1302 


EPIC-pn 
RGSl 
RGS2 


FF, 73ms 
Spectro 
Spectro 


medium 


72478 
110528 
110475 


117055 
7797 
11989 


2008 Jan 11-13 


21:32-04:07 
20:30-04:14 
20:30-04:14 


1482 


EPIC-pn 
RGSl 
RGS2 


FF, 73ms 
Spectro 
Spectro 


medium 


79549 
113617 
113846 


166180 
10022 
15921 


2009 Jan 30-31 


18:40-22:27 
18:17-22:36 
18:17-22:36 


1675 


EPIC-pn 
RGSl 
RGS2 


FF, 73ms 
Spectro 
Spectro 


medium 


72925 
101312 
101547 


185071 
11528 
17346 



<°'Net. counts in the 0.53 - 10.0 keV (EPIC-pn) and 0.35-2.1 keV (RGS) band 



May 2003 




Jan. 2008 - :f] 




Jan. 2009 




Channel energy (keV) 



Channel energy (keV) 



Fig. 1. The RGSl (black) and RGS2 (red) spectra of SN 1987 A observations in 2003, 2007, 2008 and 2009. A clear evolution in 
emission line ratios is seen. E.g. compare the nearly absent Fexvii lines at ~0.72 keV in 2003 to the strong lines seen in the later 
observations. Also the Ne x Ly„ (1.02 keV) flux increases relative to the Ne ix triplet (~0.91 keV). 



were a lready analyzed by iHaberl et al.l (l2006l) . and iHeng et aP to analyze the two new observations in combination with the 
(12008') evaluated the Jan. 2007 observation concentrating on the previous data, to obtain the evolution of the X-ray spectra as seen 
elemental abundances. In this study we use consistent models by the Reflection Grating Spectrometer (RGS. iden Herder et al] 
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I2001h . To ext end the X-ray band up to 10 keV, we included also 
the EPIC-pn ("Striider e t alJl200ll) spectra in our analysis. 

We used XMM-Newton SAS S.O.CQ to process the data. 
For the extraction of EPIC-pn spectra, single-pixel events in 
good time intervals (GTIs) with low background (threshold at 
8 cts s"' arcmin"^) were selected from the source region, cen- 
tered on SN 1987 A, and a point source free background region, 
each with a radius of 30". The RGS spectra were obtained using 
rgsproc, and GTIs (RATE < 2.0) were used to select low back- 
ground intervals. In the 2003 observation, the Honeycomb neb- 
ula is located on the dispersion axis. Using the intrinsic energy 
resolution of the RGS Focal Plane Camera CCDs, counts from 
the Honbycomb nebula are excluded. For the other three obser- 
vations, the Honeycomb nebula lies at the cross-dispersion axis 
at the very rim of the CCDs. By comparing our background spec- 
tra with the RGS background model (created by rgsbkgmodel), 
we find no significant contribution from the Honeycomb nebula 
to the background spectra. 

To obtain a detailed spectrum with ~438 ks exposure, 
the RGS spectra of the four observations were added with 
rgscombine, which also calculates combined response files. To 
allow proper application of the;t'-statistics, all spectra are binned 
to have at least 30 (RGS) or 20 (EPIC-pn) counts per bin. Due 
to the high statistics of the RGS-spectra a binning of 30 cts/bin 
does not influence most emission lines which have far more than 
30 cts in most bins, thus there is no influence on the measured 
line parameters. 

3. Spectral Analysis 

Spectral fitting was performed using XSPEC (lArnaudI [19961) 
version I2.5.0x. The errors are given for certain l^x^ ranges. 
Generally, l^x^ - 2.71 is associated with the 90% confidence 
range for one parameter of interest. 

3. 1. Identification of Emission Lines 

Since the RGS spectra of SN 1987 A are dominated by emission 
lines (see Fig.[T]), we first identified the individual lines using the 
co-added spectra. We constructed an empirical model and fitted 
it to the combined RGS 1 and RGS2 spectra simultaneously. 

As quasi-continuum we used a thermal bremsstrahlung 
model with absorption, which is formally a good approximation 
to the observed continuum, containing also radiative recombi- 
nation and two photon decays. This component also represents 
various weak emission lines, which are not resolved in the RGS 
spectra. A constant factor was allowed to vary between the two 
RGS spectra. We first fitted the strong lines with Gaussian pro- 
files. The high statistical quality of the summed spectrum re- 
quired the introduction of further weak lines where we noticed 
residua. Our final model comprises 53 lines as listed in Table |2] 
Fitted to the co-added spectra, all lines have a normalization in- 
consistent with for Ax^ = 2.71 (except for the intercombi- 
nation lines of helium-like N and Mg). For line complexes the 
line energies were combined to one value according to their ratio 
in ATOMDB 1.3.lQ and the line widths were linked. For weak 
lines it was not possible to fit individual line widths, thus we 
coupled them with the next strong line. We get a best fit with 
X^/dof = 1843/1465. The summed RGS spectra are shown in 
Fig. |2] together with the best fit model and the bremsstrahlung 
quasi-continuum (kT= 402^g eV). The line labels mark the best 

' Science Analysis Software (SAS), |http://xmni.vilspa.esa.es/sas/| 
^ http://cxc.harvard.edu/atomdb/ 
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Fig. 3. Line shifts according to the empirical model. The hori- 
zontal line shows the best fit for a systemic redshift. Errors are 
for A;^f^ = 2.71 confidence range. 



fit energy of the fitted Gaussian and name the strongest emis- 
sion line expected from plane-parallel shocked plasma emission 
codes (see below) at this energy. With these identifications we 
obtain line shifts consistent with a systemic redshift of 349 ± 24 
km s"'(see Fig.[3]l. 

The high statistical quality of the co-added RGS spectra al- 
lowed us to investigate emission lines in unprecedented detail. 
E.g. the N VI helium-like triplet at ~0.42 keV is clearly seen in 
the summed spectra (cf. Fig. [J]). Other examples of emission 
line complexes are shown in Fig. |4l Next to the Nvii Ly^ (at 
0.50 keV) line we found an additional line at slightly lower en- 
ergy (see Fig. HI, most likely N vi He^, which is the strongest 
line at this position in the plasma codes. But also argon lines 
contribute here which may cause the higher line shift. In pre- 
vious analyses with lower statistics this line likely was blended 
with the ~ 10 times more luminous N viii Ly^ line. The line shape 
at the energy of the O viii Ly^j (0.77 keV) line is inconsistent with 
a Gaussian profile and the line width for a single Gaussian is out- 
standingly higher than found for the surrounding lines (2 eV vs. 
0.25 eV). Plasma models suggest two iron lines at slightly lower 
energies, that were included in our model. The shape of the Ne ix 
triplet (~0.91 keV) is also not reproducible by three Gaussians. 
Here Fe xix (917.1 eV) may contribute to the flux, although this 
line is not expected to be strong from emission codes. 

For strong lines in the summed spectra, we find measured 
line widths inconsistent with zero within a confidence range of 
Ax^ = 2.71, demonstrating that line broadening can be detected 
significantly with RGS. Also a trend of increasing width with 
rising energy is seen. A fit of an energy dependent line width 
o-(E) = cro(£/lkeV)" yields ctq = (1.35 + 0.36) eV and a = 
2.05;^"^^. The pow er law index agrees well with the Chandra 
LEGT 2007 results (IZhekov et alJl2009l) . 

In a second step we fitted the empirical model to the spectra 
of the individual observations to follow the evolution of the indi- 
vidual emission lines in a consistent way. Only parameters with 
expected time dependence, i.e. the line fluxes and the parameters 
of the bremsstrahlung continuum, were allowed to vary. The ab- 
sorption as well as the center energy and width of the Gaussians 
were fixed at the values obtained from the summed spectra. Here 
these values can be determined more precisely than in the indi- 
vidual spectra. Using the RGS spectra with the best statistics (i.e. 
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Fig. 2. Co-added RGSl (black) and RGS2 (red) spectra from all 4 observations of SN 1987 A together with the best fit empirical 
model (green) and the absorbed bremsstrahlung continuum (blue). The line labels mark the strongest lines as expected from the two 
component vpshock model at the best fit energy of the line center. 
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Fig. 4. Examples for emission line complexes from the summed spectra: N vii Ly „, O vm Ly^j and the Ne ix triplet as seen by the 
RGSl (black) and RGS2 (red). Also the empiric model and its components are plotted. In the right panel an indication for a Fe K 
feature is seen in the summed EPIC-pn spectrum of the newer three observations. A sum of the high temperature component plane- 
parallel shock models derived from the individual observations (black) is compared with a model containing an additional Gaussian 
(red). Horizontal axes: Channel energy in keV, Vertical axes: Counts s ' keV '. 



the 2009 data), we compared line widths fixed in the individual 
fits with those allowing line widths free in the fits. We found that 
20 of 21 line widths are consistent. A time dependent analysis 
of the line widths is done in section [X21 by using plasma mod- 
els and assuming a power law dependence. The line fluxes with 
fixed and variable widths also are consistent and we conclude 
that the fixed Une widths do not influence the derived Une fluxes. 

The individual line fluxes are listed in Table |2] light curves 
for prominent oxygen and neon lines are shown in Fig. |5] 



We note, that the temperature of the quasi-bremsstrahlung con- 
tinuum shows a trend of an increase (254^^^, 413;^2i' 436;^jg 
and 405^11 for the 2003, 2007, 2008 and 2009 observation, 
respectively), but stress, that this component not represents 
bremsstrahlung only. 

We also searched for a Fe K line complex between 6.4 and 
6.7 keV in a summed EPIC-pn spectrum of the three observa- 
tions between 2007 and 2009, where we have a better statis- 
tics in this energy band. To characterise the feature, we fit- 
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Fig. 5. The lower panels show the light curves of various emission lines with fluxes in 10"^ photons cm"^ s The upper panels 
show the corresponding line ratios (not absorption corrected). The G-ratio is defined as the helium like line flux ratio of (/ + i)/r. 
Overlapping data points are shifted by 10 days for clarity. 



ted a bremsstrahlung continuum (kT=(3.3 + 0.5) keV, EM — 
(1.3 ± 0.3) X lO^'^cm-^) and one Gaussian in the 4.0-10.0 keV 
band. 

The best fit Gaussian centre energy is (6.57+0.08) keV with 
a sigma-width of 62.7+j;^i gV and the line flux is 7.3+" x 10"^ 

° -62.7 -4.2 

photons cm"^ s '. Already. iHeng et all ( l2008h noted a possible 
detection in the 2007 data, but the individual EPIC-pn and the 
summed EPIC-MOS spectra have not enough statistics for a de- 
tailed analysis. With increasing ionization, the centroid energy 
of the Fe K line complex is shifted from the 6.4 keV fluores- 
cent line over various ionization stages to the He-like emission 
of Fexxv at ~6.7. To account for this effect, we fitted an ex- 
posure weighted sum of the plane-parallel shock model derived 
for the high temperature component in Sect. 13.21 but used NEI- 
model version 1.1, which contains also Fe-ions less ionized than 
He-like. We only allowed a constant factor to account for the 
normalization and obtained a flux of the iron K lines from the 
shocked plasma of 4.9 x 10"^ photons cm"^ s"' . As can be seen 
in the right panel of Fig. |4] the resulting line shape does not 
explain the observed feature at the lower energy side. We inves- 
tigated the possibility of an additional line with zero width (also 
shown in Fig.|4|. The best fit values are 4.7 + 3.3 x 10"^ photons 
cm"^ s"' at a centre energy of 6.48;|;q }2 keV. 



3.2. Plasma Evolution 

To probe the physical parameters of the X-ray emitting plasma, 
we fitted plasm a emission models. Neithe r in XM M-Newton 
(e.g. lHeng et"al]|2008.) . nor in Chandra (e.g. lZhekov et al,.2006,) 
observations, the spectrum was described well by just one emis- 
sion component. Zhekov et al. derived a bimodal tempera- 
ture distribution in differential emission measure, peaking at 
~0.5 keV and ~2.5 keV. A non-t hermal corapone nt might also 
contribute to the hard component dPark et al.ll2005h . We noticed, 
that the emission lines in the XMM-Newton spectra can also be 
described by the sedov-model in XSPEC, which has a tempera- 
ture distribution according to the Sedov self-similar solution for 
SNRs. Here a power law is necessary for the high energy tail (cf. 
Fig.|6]). 



The abundances of the sedov component are significantly 
higher than normally found in X-ray analyses, thus we doubt the 
physical correctness of this component, but the model demon- 
strates the possible contribution of a non-thermal component up 
to (1.44 + 0.03) X 10"^ photons keV^'cm^^s"' at 1 keV with 
power law index F = 3.00 + 0.02 (for the 2009 observation). 
Formally the sedov+power law model results in a better fit (e.g. 
for the 2009 observation: = 1.35 vs. 1.48). 

However, we decided to use the approximation of a plane- 
parallel shock structure at two temperatures as commonly found 
in the literature, since the sedov model is numerically slow and 
the RGS does not provide high resolved lines above 2.0 keV. 

For a precise determination of the high temperature plasma 
component, we extended the energy band up to 10 keV by fit- 
ting the RGS 1,2 and EPIC-pn spectra simultaneously, as it is 
successfully done for calibration (Plucinskv et al. 2008). 

A comparison of the empiric model fitted to the 2009 RGS 
data and to a dataset containing also the EPIC-pn spectrum in 
the 0.2-2.0 keV energy range shows, that again 20 of 21 line 
widths were consistent. Because of the lower energy resolution 
of EPIC-pn the line widths are determined by the RGS spectra. 
We found inconsistent fluxes of the helium-like Nvi and Nvii 
Ly„ line between the EPIC-pn and RGS spectra. This is most 
likely due to calibration problems in the modeled redistribution 
of the EPIC-pn response at low energies. Thus we decided to 
limit the EPIC-pn band to 0.53-10 keV for this study. In future 
studies the 0.2 - 0.53 keV EPIC-pn band with advanced calibra- 
tion will yield additional information. 

To derive the plasma evolution we fitted a two componen t 
plane-parallel shock model (vpshock, iBorkowski et al.l 1200 ih 
with NEI-vers ion 2.0 and io n ization timescale range < t < Tu. 
Analogous to IZhekov et al.l (l2009i and references therein) the 
chemical composition of the plasma component is assumed to 
be constant in time and the abundances were fitted with the ex- 
ception of He (set to 2.57), C (0.09), Ar (0.54), Ca (0.34) and Ni 
(0.62), which produce no significant features in our data. The 
photo-electric absorption by the Galactic interstellar medium 
was set to Nh = 6xl0^"cm"^, whereas the LMC column den- 
sity with abundances set to 0.5 for metals was a free param- 
eter The abundances for absorption and e mission components 
are given according to IWilms et all (l2000h . which resulted in a 
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Table 2. Identified lines in the RGS spectea. 



Line 


E^"\ (keV) 

rest 


(keV) 

obs ^ ^ 


cr<^' (eV) 


detected flux''" (10" 


photons cm 


-2s-') 










2003 


2007 


2008 


2009 


N Vlf 


0.4198 


0.4193 


0.01 


<ii.c> 


< JO. / 


<31.4 


20.3 ±8.7 


N Vli 


0.4263 


0.4257 


0.01 


<:>4. / 




<8.7 


<13.1 


N VIr 


0.4307 


0.4301(2) 


0.01 (0-0.37) 


^1 1 A 

<:>i.4 


o /I o _i_ 1 n "7 
z4.y ± lU. / 


29.1 ±9.9 


22.8 ±9.0 


NVIHeyS 


0.4980 


0.4965(3) 


0.02 (0-0.74) 


/.o ±o.y 


y.4 ±4. J 


15.4 ±5.0 


10.6 ±5.2 


NVIILya 


0.5003 


0.4997(1) 


0.37 (0.26-0.46) 


Ib.l +l.y 


linn 1 o < 


122.4 ±9.0 


134.3 ±10.1 


Vllf 


0.5611 


0.5604 


0.40 


19.9 ±9.6 


28.9 ±5.6 


37.9 ±8.9 


51.3 ±9.8 


Vlli 


0.5686 


0.5680 


0.40 


1 n Q _i_o n 


<1 J.J 


16.1 ±10.6 


17.5 ±10.1 


O Vllr 


0.5739 


0.5733(1) 


0.40 (0.04-0.65) 


13.1 ±8.4 


43.5 ±8.8 


53.1 ±12.2 


51.3 ±10.0 


N VII Ly /? 


0.5929 


0.5926(3) 


(0-0.77) 


<Vi.i 


ij.4 ±o.Z 


22.6 ±6.4 


23.3 ±7.6 


N VII Ly y 


0.6254 


0.6247(7) 


0.63 


<y.u 


<o. / 


<5.3 


5.2 ±4.7 


N VII Ly 6 


0.6404 


0.6406(8) 


0.63 


ICQ 1 C T 

Ij.i ±J.I 


<1. / 


<2.3 


7.2 ±4.6 


O VIII Ly a 


0.6537 


0.6529(1) 


0.63 


4U.3 ±5.8 




167.5 ±8.4 


177.8 ±9.9 


OVIIHeyS 


0.6656 


0.6650(5) 


0.63 (0.46-0.72) 


y.4 ±4.y 


14. z ±4.y 


8.2 ±4.9 


16.3 ±5.5 


VII Hey 


0.6978 


0.697(2) 


0.11 


< 1 1 A C 
J. I ±4.3 


1 

<4.1 


6.9 ±4.3 


<2.8 


Fe XVm 


0.7035 


0.703(1) 


0.11 


<Z.Z 


<o.2 


<7.4 


9.8 ±4.4 


VII He 6 


0.7127 


0.712(1) 


0.11 


<8.8 


1 1 
±o.l 


5.3 ±4.0 


<3.6 


Fe XVII 


0.7252 


0.7250 


0.11 


<li. 1 


C/1 1 1 1 A T 

j4. 1 ±1U. / 


90.9 ±12.5 


84.9 ±13.7 


Fe XVII 


0.7271 


0.7269(1) 


0.11 (0-0.54) 


<o.U 


OA A 1 A /I 

zy.y ±y.4 


28.2 ±10.1 


46.1 ±11.7 


Fe XVII 


0.7389 


0.7380(2) 


0.25 (0-0.72) 


j.U ±4.4 


j4.4 ±r>.u 


38.9 ±5.2 


49.2 ±6.0 


FeXIX 


0.7696 


0.7658 


0.20 


<y.i 


CI 1/1/1 

j.i ±4.4 


<9.4 


<10.1 


Fe XVIII 


0.7715 


0.7708 


0.20 


0./ ±J.D 


iz.j ±J.Z 


24.9 ±5.9 


21.5 ±6.6 


O VinLy/? 


0.7746 


0.7738(2) 


0.20 (0-0.86) 


<y.3 


JJ.D ±J.O 


36.2 ±6.1 


50.3 ±7.1 


Fe XVm 


0.7835 


0.7817(8) 


0.25 


<y.4 


O. / ±4.1 


5.6 ±4.3 


7.7 ±5.1 


Fe XVm 


0.7935 


0.7929(9) 


0.25 


^'1 A 


1 

</.l 


8.6 ±4.0 


5.0 ±4.5 


FeXVn 


0.8023 


0.803(1) 


0.25 


<3.D 


1 1/1/1 

Iz.o ±4.4 


5.7 ±4.7 


<9.9 


Fe XVII 


0.8124 


0.8113 


0.25 


<5.2 


26.6 ±5.9 


44.9 ±6.4 


45.3 ±7.2 


O VIII Ly 7 


0.8169 


0.8159(3) 


0.25 (0-0.78) 


D.i ±4.0 


ICC 1 /C 1 

O.J ±o.l 


19.8 ±6.6 


27.3 ±7.4 


Fe XVII 


0.8258 


0.8248(1) 


0.64 (0.23-0.89) 


14. y ±4. J 


5Z. / ±o.4 


100.3 ±6.9 


123.5 ±8.0 


O VIII Ly 6 


0.8365 


0.835(2) 


1.89 


<4.y 


^ A C 

<4.j 


5.2 ±5.0 


12.5 ±5.6 


Fe XVIII 


0.8530 


0.8517(8) 


1.89 


<i.L 


11-1 1 /I 

1 1.4 ±4.0 


16.5 ±5.0 


19.4 ±5.7 


Fe XVIII 


0.8626 


0.862(1) 


1.89 


<1U.U 


1 1 Q _i_C c 

11. J ±J.J 


18.3 ±5.5 


13.4 ±6.0 


Fe XVIII 


0.8726 


0.8719(4) 


1.89 (1.45-2.35) 


<y.o 


QA T 1 C 1 
jU. / ±J.l 


39.9 ±5.5 


54.6 ±6.5 


FeXVn 


0.8968 


0.8957 


0.85 


D.O ±4./ 


MA 1 /I C 

lj.4 ±4. J 


17.1 ±4.9 


14.5 ±5.5 


NeEXf 


0.9050 


0.9039 


0.85 


17.0 ±6.6 


/IOC 1 T A 

43.5 ±7.0 


61.8 ±7.9 


57.7 ±8.8 


NeEXi 


0.9148 


0.9136 


0.85 


<il.O 


1/1 _i_n 1 
14. o ±12.1 


23.2 ±13.7 


<30.1 


FeXIX 


0.9171 


0.9161 


0.85 


<13.2 


<2j. / 


<26.3 


38.2 ±19.0 


NelXr 


0.9220 


0.9208(2) 


0.85 (0.25-1.23) 


ion 1 T o 

12. y ± /.6 


O C Q 1 1 1 A 
6J. J ± 1 l.U 


92.7 ±11.3 


99.7 ±12.8 


FeXX 


0.9651 


0.965(1) 


1.76 


^A Q 

<4. J 


1 A A _i_C O 

14.0 ±J.o 


20.1 ±6.2 


24.1 ±6.9 


Fe XVII 


1.0108 


1.009(1) 


1.76 


^n o 

<y.8 


1 /.o ±o. / 


28.1 ±7.6 


17.7 ±8.0 


Ne X Ly a 


1.0219 


1.0205(3) 


1.76(1.38-2.28) 


zU.8 ±/.i 


nc 1 _i_n o 

yj.j ±y.2 


107.0 ±9.9 


145.3 ±11.5 


Ne IX He yS 


1.0740 


1.075(2) 


3.76 (2.41-5.53) 


n A _i_o n 
y.4 ±8.U 


1 C _i_'7 ^ 
lO.J ± /.J 


24.9 ±7.3 


24.5 ±8.0 


Ne IX He 7 


1.1270 


1.124(3) 


1.34 


A 

<i A 


1 1 O 1 C 1 

1 1.2 ±j.l 


10.5 ±5.7 


7.0 ±5.9 


FeXVn 


1.1512 


1.15(5) 


1.34 


<j.l 


y.O ±J.O 


10.4 ±4.3 


13.3 ±4.7 


NeXLyyS 


1.2109 


1.210(2) 


1.34(0-3.94) 


<10.9 


<9.4 


<8.3 


<8.8 


MgXIf 


1.3311 


1.3297 


2.25 


9.0 ±8.3 


10.7 ±5.9 


13.9 ±5.7 


18.6 ±5.9 




1 ^zL^ 1 


1 'XAM 
1 .^H-l / 




<11.8 


<13.0 


<.o.o 


< / .u 


Mg XI r 


1.3522 


1.3508(9) 


2.25 (0-3.73) 


<15.3 


25.9 ±7.0 


35.2 ±6.7 


44.4 ±6.6 


Mg XII Ly a 


1.4726 


1.473(2) 


4.25 (0-7.27) 


<15.0 


12.7 ±4.5 


21.1 ±5.4 


26.6 ±5.7 


Mg XI He yS 


1.5793 


1.582(6) 


5.01 (0-17.1 ) 


12.6 ±10.0 


5.6 ±4.5 


7.4 ±5.3 


6.2 ±5.1 


Si XIII f 


1.8394 


1.840 


5.38 


<31.0 


<20.9 


<27.4 


<27.9 


Si xm i 


1.8537 


1.855 


5.38 


<33.9 


<29.3 


<4.96 


<38.7 


Si xm r 


1.8649 


1.865(4) 


5.38 (0-15.8 ) 


<36.3 


34.4 ±11.1 


<4.68 


39.8 ±12.6 


Si XrV Ly a 


2.0060 


1.997(1) 


(0-230) 


<85.6 


22.5 ±12.4 


19.7 ±19.4 


18.1 ±15.8 



Note: All errors are for Ax^ = 2.71 confidence range. Parameters without errors were coupled to the corresponding parameter of the following line 
with higher energy (sec text). 

Rest energy of identified line according to ATOMDB 1.3.1. If the line is a multiplet, we give the energy of the strongest component. 
**' Central energy of the fitted Gaussian. 

Best fit value of Gaussian sigma width. 
'* Detected flux for the individual observations. For fluxes consistent with at Ax^ < 2.1 \ only the upper Umit is given. 
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Fig. 6. The XMM-Newton RGSl (black), RGS2 (red) and EPIC-pn (green) spectra of the 2009 observation. Plotted is also a model 
containing two plane-parallel shock components (left panel) and a sedov+powerlaw model (right panel). 



Table 3. Results from the VPSHOCKh-VPSHOCK model fits. 



Parameter May 2003 Jan. 2007 Jan. 2008 Jan. 2009 



/tri(keV) 


41+0-'" 

"•^^-0.01 


r\ CO +0.002 
"•-^ -0.003 


n CQ+0.002 
"•-^-'-0.001 


n C4+0.002 
"•-^ -0.001 




2 89+° " 


7 4^+0.06 


2 50+°-°^ 


a o +0.04 




4.09:°;! 


17.83!°:!^ 


23 45+° '^ 


28.37!°;!^ 


EMf 


1 06+°"^ 


7 qc+0.04 
^■'^■^-o.os 


T q 1+0.06 

-^•^^-0.05 


4 q 1+0.05 
^•^^--0.07 


u,l 


T 9^+0.29 




6 03+° ''' 

"•"-'-0.14 


6.86^?^ 


_(« 

u,2 


1 TO+0.12 


1 q^+o.ii 


2.02!°:?" 


9 99+0.13 
^•^^-0.09 


o-o(eV) 


3 51+'-" 


1 51+° 25 


1 SS+°-31 
-^"-0.29 


1 43+°" 
-^-0.29 


a 


1 10+°-^'' 

^•^"-0.76 


2 64+° ® 

-0.45 


T 1 c+0.43 
-'■^•'-0.63 


3 10+°'' 
-'•^"-0.52 



Note: Errors for = 2.71 confidence range. 

" EM= J ngfiHclV in 10^* cm"', a distance of 50 kpc is assumed. 

* T^: upper limit on the ionization time range in 10" cm"' s. 

slightly better fit (y ^ , - 1 .32 vs. 1.38) than the abundances of 
lAnders & Grevessa ( 19891) . The main impact is on the absorp- 
tion, which influences continuum emission and line ratios. We 
found no impact on the main results of our study. A constant 
factor is allowed to vary for calibration differences of the indi- 
vidual instruments. We modified the shock model line widths to 
be described by the power law function mentioned above, thus 
our custom version of the vpshock model contains two more 
parameters (cro, a). 

We fitted the fourXMM observations (12 spectra) simultane- 
ously, which results in a;i^^/dof = 5400/4084. As best fit param- 
eters we obtained for the LMC absorption Nh= (2.23 + 0.02) x 
10^' cm"^, for the systemic velocity vgvA - 326![g km s"'and 
for the instrument dependent constants Crgs2 = 1.03 + 0.01 and 
Cpn = 1.04 + 0.01 (relative to Crgsi = !)■ For the time depen- 
dent parameters see Table [3] and Fig. |7] for the abundances see 
Tables 

3.3. Integrated Flux 

To derive the dete cted flux analogous to lHaberl et al.l ( l2006h and 
iHeng et aP (l2008h . the vpshock+ vpshock model was fitted sepa- 



rable 4. Abundances of SN 1987 A 



2vpshock sedov+po HOB Z09 LF96 H98 



N 


1 ^cs+°°'*'' 


o -34+0.34 


0.54 


0.83 


2.37 







n 1 98 +0.002 
"■^^°-0.002 


40+° °' 
"■^ -0.01 


0.05 


0.14 


0.33 


0.33 


Ne 




1 00+°°' 


0.20 


0.41 


0.63 


0.41 


Mg 


n 9Q 1 +0.008 
"•^'^-0.007 


7J+0.04 

' ^-0.03 


0.17 


0.42 




0.48 


Si 


516+°°'** 


1 01+°°** 

^■"^-0.04 


0.48 


0.63 


0.91 


0.59 


S 


45 J +0.039 
"•^-'^-0.041 


1 56+°-' 


0.59 


0.40 


0.46 


0.48 


Fe 


n 994+0.003 
" -0.002 


46+°'° 
"■^"-0.01 


0.09 


0.33 


1.12 


0.38 



Note: Abundances relative to IWilms et all (l2000l) . 2vpshock: to com- 
ponent plane parallel shock model fitted simultaneously to all 12 spec- 
tra (this study), sedov+po: sedov model with non thermal component 
fitted to the 2009 data (this study), H08: EPIC-pn 2007 dHeng et all 
l2008h VPSHOCK+ VPSHOCK WOO model, H09: Chandra result s 
dZhekov et al ] |2009l). LF96: Inner ring dLundgvist & Fransson|[T99^ . 
H98: LMC-average dHughes et all 19981) . 

rately to the most recent two EPIC-pn spectra. The fluxes do not 
depend strongly on individual model parameters. The fluxes for 
various sub-bands are given in Table|5] 

4. Discussion 

We monitored SN1987A with XMM-Newton and obtained 
RGS spectra with unprecedented statistical quality. We analyzed 
individual line shapes and fluxes with an empirical model and 
applied physical models to a combined set of RGS and EPIC-pn 
spectra. 

The line shifts derived from the empirical model (349 ± 24 
km s"') and the plasma codes (326 + 16 km s"') are con- 
sistent. Considering the RGS wavelength scale calibration the 
line shifts are consistent with the systemic velocity derived 
from optical SN 1987 A observations (286.74+0.05 km s"', e.g. 
iGroningsson et al.ll2008l) . 

Line broadening is significantly detected in the RGS spec- 
tra. To obtain a conversion factor for the measured line widths 
and the bulk gas velocity, we convolved the emission spectrum 
of a cylindrically expanding ring at an inclination of 45° and 
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7. Time evolution of the plasma variables: Horizontal axes are days after explosion. T] 2 is given in keV, EMi 
in 10"cm"^s and ctq in eV. 



8000 

,2 in lO^'^cm"^ 



Table 5. EPIC-pn fluxes 



sub-band 
(keV) 



Rev. 1482 
Jan 2008 



Rev. 1675 
Jan 2009 



0.2- 


-0.8 


1 QO+0.01 


1 sq+o.oi 
^■•^'-0.01 


0.8- 


-1.2 


2 04+°"- 


2 48+°"' 

^•^°-0.02 


1.2- 


-8.0 


2 06+°"^ 


9 cQ+0.05 
-^■-^-0.04 


0.5- 


-2.0 


4 -i 1 +0.03 


=; 94+0.04 

-^•^ -0.04 


3.0- 


10.0 


51+°"' 
"■-^ -0.02 


n ss+0.02 


0.5- 


10.0 


5 90+0.05 

°-0.08 


6 40+°"^ 


0.2- 


10.0 


c CI +0.06 
•^■•^^-0.07 


6 64+°"' 



Note: Fluxes in 10 '^ erg cm ' s '. Errors are f or 68% confidence. 
For the fluxes of t he previous observations see iHaberl et al.l d2006) and 
iHeng etalj | |2008h . 



a temperature of kT - 0.5 keV in a bulk gas velocity range 
of < v'buik < 1000 km s ' with the RGS response and fit- 
ted a Gaussian to the resultant line profile. In this simulation 
we get Vbuik = 0.80 x Av (FWHM). For a similar analysis see 
iMichael et al](l2002h . 

As seen in Fig. [8] Vbuik is determined best by the profiles of 
the two strongest lines N vii Ly„ and O viii Ly„. The broaden- 
ing of the former line including a statistical error of Ax^ = 2.71 
is SlS+Jjy km s"' FWHM. Since the instrumental line broaden- 
ing is somewhat higher, the uncertainty of the RGS line spread 
function (Isf) causes a similar error A 10% uncertainty of the Isf 
(A/lisf ~ 0.07A FWHM) gives an uncertainty range of 400-600 
km FWHM. But recent studies indicate that the current Isf- 
model rather over-estimates the instrumental line width (Jelle 
Kaastra, private communication). For our modeling of the RGS 
lines the extent of the source (~1'.'5) is negligible. For a ther- 
mally equilibrated plasma, as it is expected fo r the shock veloc - 
ity range of the lower temperature component (lRakowskill2005h . 
the thermal line broadening also does not contribute significantly 
(e.g. 133 km s^' FWHM for nitrogen at kT = 0.5 keV). We also 
found no influence due to the two unresolved Lyman transitions 
by modeling the line with two Gaussians, and other unresolved 
lines should not contribute significantly either 

For the Nvii Lyo- line this indicates a bulk velocity of 
(414 + 120) km s"'. In the case of gas shocked by a strong adi- 
abatic shock wave and thermal equilibration, this yields Vsh - 
4/3vbuik - 552 km s^'and a post shock temperature of kT - 



co-added 
May 2003 
Jan. 2007 
Jan. 2008 
Jan. 2009 




Energy (keV) 

Fig. 8. Squares denote the bulk gas velocity inferred from the 
emission line broadening derived from the Gaussian model fit- 
ted to the co-added RGS spectra. Errors are for Ax^ = 1 . The 
solid black line shows the best fit power law function to these 
data points. The other (red) lines show the best fit power law 
functions as derived with the modified vpshock model from the 
individual spectra. For clarity, the errors are not plotted. Note, 
that especially for the 2003 data the eiTors are large (cf. Fig. Q 
right panels). The horizontal blue lines indicate the post shock 
gas velocities resultant from strong adiabatic shocks that heat 
the plasma to a temperatures of kT - 0.5 and 2.5 keV after ther- 
mal equilibration. 



3/16K 



0.435+n ?T keV. Here an adiabatic index of - 



•sh - "■-^-^•'-0.22 "^"^ '■ "" — — r - 5/3 
and an average molecular weight of /v = 0.73mp according to the 
abundances of the inner ring is assumed. Similarly for the O viii 
Lyo- line we obtain kT = 0.75+^ 3^ keV. Thus these line widths 
are consistent with the temperatures derived from the plasma 
models. 

The bulk gas velocities derived from these two lines are also 
consistent with the result of the Chand ra LETG observation in 
2007 (-360 km S-' , IZhekov et alj2009h . But in general the bulk 
gas velocities as function of line energy which can be described 
globaly by power laws, are higher for the RGS derived values 
than the values deduced from the Chandra spectra (e.g. by a fac- 
tor of ~2 at 1 keV). Some of this difference is caused by thermal 
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broadening, which contribute to the XMM-Newton derived line 
widths, but not to the Chandra values, which are based on spatial 
spectral deformation (i.e. the bulk gas velocity is deduced from 
the difference of the line broadening in the two grating arms). 
But this effect would only become significant for plasma out of 
thermal equilibrium. Also the radial velocity distribution of the 
shocked plasma might contribute differently to the line broaden- 
ing seen in the RGS and dispersed Chandra spectra, but it can 
not explain such a large difference. We believe that most of the 
difference is caused by the different modeling methods and by 
calibration differences. 

With our modified plasma emission code we can also inves- 
tigate the evolution of the line broadening and we can narrow 
the artificial broadening due to blended lines. The increase in 
a indicates that the velocity of regions with lower temperature 
decreases relative to the regions with higher temperature. The 
RGS derived line widths are determined most accurately in the 
0.5-1.0 keV band, where we have the best statistics and strong 
emission lines. Conspicuous is the sudden decrease of the line 
broadening in this energy band between 2003 and 2007, indi- 
cating a deceleration of the bulk velocity. Although the statistics 
of this observation is low, a comparison with a line broadening 
fixed to the function derived from the empirical model shows 
that the individually derived line description is significant (f-test 
probability 1.4 x lO^^). From 2007 to 2009 the XMM-Newton 
derived line widths show a rather constant line width, similar 
to those of the Chandra grating observati ons in 2004 and 2007 . 
But from the HETG observation in 1999 iMichael et all d2002h 
report line widths corresponding to a post shock plasma velocity 
of -2500 km s Thus the line broadening had to decrease in 
between the XMM-Newton observation (5918 SN days) and the 
Chandra LETG sequence (-6400 SN days) by -50% at 1 keV. 
We note, that the sudden deceleratio n of the ex pansion seen in 
the X-ray images around day 6100 dRacusin et al. 2009) also 
falls into this time interval. But we emphasize that the spectra 
in 2003 might be more complex than assumed by two shock 
components and could contain a more complex mixture of ra- 
dial velocity and temperature components. 

Also the derived plasma parameters of the plasma emission 
model show a clear evolution between 2003 and 2007. The emis- 
sion measure ratio EM2/EM1 shows a decrease. During this 
time also the soft flux was rising drastically. In the later obser- 
vations the EM ratio is rather constant with a tendency of an 
increase. 

We also observe the ongoing ionization of the plasma. This 
causes e.g. the increase of the [Ne x Ly„]/[Neix He r] flux ratio 
by -26% during the last two observations. The upper limit of 
the ionization timescale range Tu is steadily increasing for both 
plasma components. Tu 1 indicates some upturn after the first ob- 
servation, which could be caused by an increasing density. A fit 
of a linear function to the last three values (see Fig.|7]) yields an 
"average" density of rii = (3.24 ± 0.64)xl0"' cm"-'. For the high 
temperature component no evolution in density is observed and 
a similar fit yields «2 - (4.8 + 1.4)xl0^ cm"-'. The ratio of the 
emission measure of the low and high temperature components 
is roughly constant with a volume ratio of -7.5. 

The densities derived from the unshocked gas in the optical 
are 6 x 10-* cm"^ < n < 3 x lO'* cm"^ for the ring and - 10^ 
cm"^ for the extended nebula dLundqvist & Franssonlll996h . A 
strong adiabatic shock wave increases these densities by a fac- 
tor of four Our values would suggest less denser regions in the 
ring and the surrounding nebula as origin of the X-rays. But the 
vpshock model assumes a ionization time distribution linear in 
emission measure, which is clearly not true for SN 1987 A. A 



complex density distribution and effects of radial expansion af- 
fect the derived densities. However, these values are not consis- 
tent with those from the m odeling of the light curve (Haberl et al] 
l200a lAschenbachll2007h . and should be regarded with some 
care. 

The temperature of the soft plasma component shows a 
strong increase between the first two observations, and a slight 
increase in the later ones. Also from Chandra observations a ris - 
ing temperature of this com ponent is de r ived (Park et al.ll2006l) . 
Between 2004 and 2007 .Zhekov et all (120091) report a nearly 
constant temperature. With RGS we can see an ongoing small 
increase in the temperature of the soft component. However, a 
decrease caused by the deceleration of the shock wave and adia- 
batic expansion of the shocked plasma is expected which is seen 
in the high temperature component. An increase is possible, if 
regions with slightly higher temperature contribute more to the 
emission measure with time, e.g. due to more rapidly rising vol- 
ume of these regions. But also the effects of incomplete thermal 
equilibration might contribute to a rising electron temperature, 
especially initially after the shock wave has reached denser pro- 
trusions of the ring. 

The abundances derived from the two component plan- 
parallel shock model (cf. Table HJi are higher (on aver- 
age by a factor of 1.8), than derived by iHeng et al] (l2008l 
VPSHOCK+ VPSHOCK WOO model). Heng et al. used a sim- 
ilar model, but fitted only the EPIC-pn spectum of the 2007 
observation and derived different plasma parameters (e.g. tem- 
peratures) which influences the derived ab undance s. Com pared 
with the Chandra grating results of Zhe kov et al. 1 (l2009h . our 
abundances are consistent within 30%, except nitrogen, where 
our value is higher (by a factor of 1.7) but rather matches the 
iLundqvist & FranssonI (1 1996) value, and Mg and Fe, which are 
-30% lower. We see, that the individual modelling causes sys- 
tematic differences, much larger than the statistical errors. We 
also note, that the abundances derived from the plane-parallel 
shock model are lower than the abundance s of the inner ring as 
derived from ILundqvist & FranssonI dl996t) . whereas the abun- 
dances of the sedov-Hpowerlaw model are higher Thus the lower 
abundances derived so far in X-ray analysis may be caused by 
the assumption of a two component plan-parallel shock struc- 
ture. 

As shown in the right panel of Fig. |4] we see a clear in- 
dication for an Fe K feature. The black curve in Fig. |4] shows 
the emission from the shocked plasma of the inner ring, as ex- 
pected from our two component shock model, which describes 
the Fe L lines well. Surprisingly we find, that the feature sug- 
gests a contribution of emission at lower energy of less ionized, 
or even neutral, iron as it is present in the unshocked part of the 
ring and in the supernova debris. In the case of a very young 
(Tu < 3 X 10'") reverse shock in a iron rich region, we would 
also expect emission from other elements in the RGS spectra. If 
the additional emission is caused by fluorescence (6.4 keV), this 
would need neutral iron. In this case, the iron could be excited 
by X-rays from the shocked ring, but this would need a much 
higher column density of the unshocked region. Thus the emis- 
sion might originate from the supernova debris, where a higher 
column density is possible and reprocessed radiation from nu- 
clear decays might contribute. 

5. Conclusions 

1 . With our monitoring we can follow the detailed evolution of 
the supernova remnant of SN 1987 A. E.g. the upturn in ion- 
ization age and emission measure ratio between 2003 and 
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2007 shows that the blast wave was propagating into the in- 
ner ring. 

2. The decreasing Une widths at lower energies in between the 
first two observations indicate a deceleration of the lower 
temperature plasma, that correlates well with the decelerat- 
ing ring expansion as observed in the Chandra images. 

3. The electron temperature derived for the soft temperature 
component with plasma models is consistent with the line 
widths of the emission lines in the corresponding energy 
range. This is expected, if the emission is primary caused 
by shocks transmitted into denser regions. 

4. The lower statistics at higher energies do not justify such 
conclusions for the high temperature component, where the 
bulk velocity can be reduced by the contribution of reflected 
shocks which would also heat the plasma further. But the 
decreasing temperature with time and the higher Une widths 
seen by XMM-Newton and Chandra rather suggest forward 
shocks in less denser regions as the dominating process. 

5. The iron K feature argues rather for a thermal high energy 
component and little if any non-thermal contribution. The 
line shape further suggests the contribution (~50%) of a cold 
iron line (23cr level), possibly emitted from the supernova 
debris. 

Further monitoring will yield information on the arising SNR 
of SN1987A and the increasing statistics of the brightening 
source will allow even more detailed analyses also beyond the 
assumption of plane-parallel geometry. This might put further 
constraints on the origin of the X-ray emitting regions and their 
dynamics. Similarly the iron K emission can yield information 
on the evolution of the debris. 

References 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, 

Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. 

Jacoby & J. Barnes, 17 
Aschenbach, B. 2007, in American Institute of Physics Conference Series, 

Vol. 937, Supernova 1987A: 20 Years After: Supemovae and Gamma-Ray 

Bursters, ed. S. Immler, K. Weiler, & R. McCray, 33^2 
Beuermann, K., Brandt, S., & Pietsch, W. 1994, A&A, 281, L45 
Borkowski, K. J., Lyerly, W. J., & Reynolds, S. R 2001, ApJ, 548, 820 
den Herder, J. W., Brinkman, A. C, Kahn, S. M., et al. 2001, A&A, 365, L7 
Dewey, D., Zhekov, S. A., McCray, R., & Canizares, C. R. 2008, ApJ, 676, L131 
Gaensler, B. M., Staveley-Smith, L., Manchester, R. N., et al. 2007, in American 

Institute of Physics Conference Series, Vol. 937, Supernova 1987A: 20 Years 

After: Supemovae and Gamma-Ray Bursters, ed. S. Immler, K. Weiler, & 

R. McCray, 86-95 

Groningsson, P., Fransson, C, Leibundgut, B., et al. 2008, A&A, 492, 481 

Groningsson, P., Fransson, C, Lundqvist, P., et al. 2006, A&A, 456, 581 

Haberl, F, Geppert, U., Aschenbach, B., & Hasinger, G. 2006, A&A, 460, 811 

Hasinger, G., Aschenbach, B., & Truemper, J. 1996, A&A, 312, L9 

Heng, K., Haberl, F., Aschenbach, B., & Hasinger, G. 2008, ApJ, 676, 361 

Heng, K., McCray, R., Zhekov, S. A., et al. 2006, ApJ, 644, 959 

Hughes, J. P, Hayashi, 1., & Koyama, K. 1998, ApJ, 505, 732 

Jansen, F, Lumb, D., Altieri, B., et al. 2001, A&A, 365, LI 

Lawrence, S. S., Sugerman, B. E., Bouchet, P, et al. 2000, ApJ, 537, L123 

Lundqvist, P & Fransson, C. 1996, ApJ, 464, 924 

McCray, R. 2007, in American Institute of Physics Conference Series, Vol. 937, 
Supernova 1987A: 20 Years After: Supemovae and Gamma-Ray Bursters, ed. 
S. Immler, K. Weiler, & R. McCray. 3-14 
Michael, E., McCray, R., Chevalier, R., et al. 2003, ApJ, 593, 809 
Michael, E., Zhekov, S., McCray, R., et al. 2002, ApJ, 574, 166 
Ng, C.-Y, Gaensler, B. M., Staveley-Smith, L., et al. 2008, ApJ, 684, 481 
Park, S., Zhekov, S. A., Burrows, D. N., et al. 2006, ApJ, 646, 1001 
Park, S., Zhekov, S. A., Burrows, D. N., & McCray, R. 2005, ApJ, 634, L73 
Plucinsky, P. P., Haberl, F., Dewey, D., et al. 2008, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7011, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Racusin, J. L., Park, S., Zhekov, S., et al. 2009, ApJ, 703, 1752 



Rakowski, C. E. 2005, Advances in Space Research, 35, 1017 

Struder, L., Brief U., Dennerl, K., et al. 2001, A&A, 365, L18 

Sturm, R., Haberl, F, Hasinger, G., Kenzaki, K., & Itoh, M. 2009, PASJ, 61, 895 

Turtle, A. J., Campbell-Wilson, D., Manchester, R. N., Staveley-Smith, L., & 

Kesteven, M. J. 1990, lAU Circ, 5086, 2 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 
Zhekov, S. A., McCray, R., Borkowski, K. J., Burrows, D. N., & Park, S. 2005, 

ApJ, 628, L127 

Zhekov, S. A., McCray, R., Borkowski, K. J., Burrows, D. N., & Park, S. 2006, 
ApJ, 645, 293 

Zhekov, S. A., McCray, R., Dewey, D., et al. 2009, ApJ, 692, 1190 



